arXiv:1505.06870v2 [math-ph] 5 Oct 2015 


ZMP-HH/15-14 

ITP-UU-15/06 

Differential approach 
to on-shell scalar products 
in six-vertex models 


W. Galleas“ and J. Lamers^ 

“ II. Institut fiir Theoretische Physik 
Universitdt Hamburg, Lumper Chaussee 149 
22761 Hamburg, Germany 

Wellington.galleasOdesy.de 

^ Institute for Theoretieal Physics, 

Center for Extreme Matter and Emergent Phenomena, Utrecht University 
Leuvenlaan 4, 3584 CE Utrecht, the Netherlands 

j.lamers@uu.nl 


Abstract 

In this work we obtain hierarchies of partial differential equations de¬ 
scribing on-shell scalar products for two types of six-vertex models. 
More precisely, six-vertex models with two different diagonal boundary 
conditions are considered: the case with boundary twists and the case 
with open boundary conditions. Solutions and properties of our partial 
differential equations are also discussed. 
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1 Introduction 

The behavior of physical systems at criticality is endowed with remarkable properties 
and its description is intimately associated with the study of correlation functions. When 
a system is away from a critical point, two-point correlation functions are expected to 
decay exponentially as the distance between the target points becomes inhnitely large. 
This asymptotic behavior motivates the dehnition of the correlation length which con¬ 
tains fundamental information concerning the system critical behavior. For instance, 
correlation lengths are expected to diverge in second-order phase transitions when a 
system reaches its critical point. In fact, this singular behavior can be regarded as a 
trademark of a second-order phase transition and this divergence is then governed by a 
power law characterized by a critical exponent. 

The scenario at criticality is drastically different and correlation functions are ex¬ 
pected to decay according to a power law instead of exponentially. This feature mo¬ 
tivates one to dehne the anomalous dimension for the corresponding order parameter. 
The anomalous dimension is also a critical exponent and this particular power-law be¬ 
havior shows that statistical fluctuations of the order parameter are strongly correlated 
throughout the entire extension of the physical system. In the vicinity of a phase transi¬ 
tion correlation lengths are much larger than the system’s lattice spacing and it can be 
regarded as a natural distance scale. Furthermore, the correlation length is a function of 
the model’s coupling constants and the configurations of the system over the correlation 
length can be made sufficiently smooth such that the continuum formalism of quantum 
field theory is appropriate. For a detailed discussion we refer the reader to [T]^. 
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The computation of correlation functions of interacting systems is a highly nontrivial 
task and it is often performed by means of perturbative expansions. However, since 
critical behavior is associated with divergences of the aforementioned quantities, it is 
highly desirable having a non-perturbative description of correlation functions. In this 
way, exactly solvable models of statistical mechanics and quantum field theory are natural 
candidates for the investigation of critical phenomena. 

In the context of quantum field theory correlation functions are usually described 
by means of differential equations. We have for instance the Callan-Symanzik equation 
arising from the renormalization group framework [4]-[^ and the Knizhnik-Zamolodchikov 
(KZ) equation describing correlation functions of primary helds in conformal held theories 
[^. Differential equations also seem to play a prominent role in the study of correlation 
functions of exactly solvable models of statistical mechanics, although results in this 
direction are still limited to a few models. For instance, in the authors have shown 
that a two-point function of the impenetrable Bose gas can be described by a differential 
equation of Painleve type. Prior to that, Painleve equations were also found in the study 
of certain two-point functions of the two-dimensional Ising model 

the 


Within the framework of the Quantum Inverse Scattering Method (QISM) 14 


study of form factors and correlation functions is intimately related to the evaluation of 
scalar products of Bethe vectors 15 . In particular, as far as the six-vertex model and the 


associated Heisenberg spin chain are concerned, the form factor of the particle-number 
operator has been written in terms of on-shell scalar products of Bethe vectors in 16 . It 


is worth remarking that the relation between form factors and on-shell scalar products 
obtained in 16 builds on a detailed analysis of form factors previously investigated 


m 


17,18 


The above discussion raises some questions concerning the description of correlation 
functions through differential equations. For instance, there is a fundamental difference 
between the differential equations satisfied by correlation functions in conformal field 
theories and those obtained for the Bose gas and the Ising model. The latter are non¬ 
linear equations of Painleve type whereas KZ equations are linear. Hence one may 
wonder if there exists linear partial differential equations (PDFs) describing six-vertex 
model correlation functions. This question is the main motivation for the present paper. 
Our arguments here shall rely on the results presented in expressing form factors of 
the six-vertex model in terms of on-shell scalar products of Bethe vectors. In this way 
we achieve a linear differential description of six-vertex model form factors by presenting 
linear PDFs determining on-shell scalar products. 

The possibility of deriving linear PDFs for such quantities was first put forward 


19 based on a detailed analysis of certain functional equations originated from the 


m 

Yang-Baxter algebra. Several quantities and lattice systems have been tackled through 
this algebraic-functional method and we refer the reader to 20,21,19,22 26 for an ac¬ 
count of the results so far. In particular, linear PDFs describing partition functions with 
domain-wall boundaries and spectral problems for transfer matrices have been presented 
in 19,24,25 and 27 respectively. Although the Yang-Baxter algebra is a fundamental 
ingredient, this approach is not limited to this particular algebraic structure. For ex¬ 
ample, integrable systems with open boundary conditions can also be included in this 
framework by replacing the Yang-Baxter algebra by the reflection algebra 


25,26 
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Outline. In the present work we construct a differential description of six-vertex models 
on-shell scalar products for two types of boundary conditions. This paper is organized 
as follows. Section is concerned with the case of diagonal boundary conditons, also 
referred to as quasi-periodic boundaries. Our analysis in this case will rely heavily on 
results previously presented in 23 which includes a functional equation for the relevant 


scalar product. These results are recalled in Section 2.1 and the corresponding PDEs 
are obtained and studied in Section |2.2[ On-shell scalar products for diagonal open 


boundary conditions are discussed in Section Here we make use of the results recently 
obtained in whose details are given in Section 3.1 
in Section 13.21 and we conclude with hnal remarks in Section 01 


The resulting PDEs are described 


2 Twisted boundary conditions 

The study of integrability-preserving boundary conditions is an important branch of the 
theory of exactly solvable models [^. For instance, it is well known that particular 
boundary terms can be used to relate the critical behavior of several models of statistical 
mechanics [1,28,29 . Relations between different systems away from the criticality are 
more subtle but nevertheless such relations still exist. This is the case for the six-vertex 
models with domain-wall and anti-periodic boundary conditions 30 , whose relation has 


been established through the same algebraic-functional approach employed in 23 . As 


far as the six-vertex model is concerned, boundary conditions are not a mere detail that 
can be disregarded in the thermodynamical limit. In that case boundary conditions 


even influence the model bulk free-energy 31-35 and modify the spectrum of finite-size 


corrections 36,29 . The latter provides fundamental data concerning the underlying con¬ 


formal held theory. We also point out that twisted boundary conditions arise naturally 
in the context of Bethe/gauge correspondence by including certain topological terms in 
the lagrangian of the two-dimensional M = (2,2) gauge theory 


37 


As far as vertex models with twisted boundary conditions are concerned, the algebraic 
formulation of the integrability-preserving case was put forward in 


38 


based on the 

invariance of the Yang-Baxter algebra under a particular linear map. This algebraic 
setting was then used in [^ to construct a functional equation governing on-shell scalar 
products. We shall build our analysis on those results in order to derive a hierarchy of 
PDEs describing such scalar products. 


2.1 Functional equation 

The starting point of our analysis is the study of six-vertex model scalar products by 


means of functional equations as described in 23 . Throughout this work we shall use 


the terminology on-shell scalar product to refer to a scalar product of Bethe vectors 
when only one of the vectors is subjected to Bethe ansatz equations. We consider these 
on-shell Bethe roots as fixed complex parameters, so that the on-shell scalar products 
can be regarded as a function iS„: C” —)■ C where n G Z>o is the number of magnons 
contained in the Bethe vectors. In its turn, the integer index n is bounded by the size 
of the rectangular lattice on which our six-vertex model is defined. In fact, if L G Z>o 
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is the lattice horizontal length, then 0 < n < L. The on-shell scalar product Sn is a 
multivariate function in n complex variables referred to as spectral parameters. Before 
proceeding let us pause for a moment to set up our notation. 

Definition 1 (Indices, variables and parameters). Let n G Z>o satisfy the condition 
0 < n < L for L G Z>o. Then consider indices k E Z and let A^, Af G C denote 
spectral variables. In addition to that, we shall use 7 G C to denote the anisotropy 
parameter, G C will be inhomogeneity parameters and (fi, 02 £ C\{0} will characterize 
the boundary twists. For convenience, we shall also employ the short-hand notation: 
a(A) := sinh(A -|- 7 ), 6 (A) := sinh(A) and c(A) := sinh( 7 ). 

Definition 2 (Spectral parameter and inhomogeneity sets). Define A*’-^ := {A^ G C | z < 
k < j} and := A®’-^\{A}. It is also convenient to define the inhomogeneity set 
< k < L} in order to recast the dependence of the six-vertex model 
scalar products with the inhomogeneity parameters. 

Definition 3 (Bethe roots). Let Af G C for 1 < k < n. We then define the set 
B := {Af, Af,..., A^} of spectral variables solving 

~ a) ^ TT ~ 

K^k - l) (1>1 fj- a(Af - Af) ■ 

i^k 

The set B will also be referred to as Bethe roots set. 

Remark 1. The constraint imposed on the variables Af, i.e. Bethe ansatz eguations, 
admits several inequivalent sets of Bethe roots. Througout this work we let B be any 
single fixed set of such Bethe roots. 

Among the main results of we have the following theorem. 

Theorem 1 (Functional equation). Taking into account the above definitions, the on- 
shell scalar product Sn satisfies the functional equation 

Ko Sn{A^n + 5„(A°/) = 0 , ( 2 . 2 ) 

AeAi." 


n 


with coefficients Kq and K\ given by 


Ko := 01 JJ a(Ao - /z) 


fi^Vl 


n 


a(A — Ao) 


n 


a(A — Ao) 


, KA-Ao) fi6(A-Ao) 

LasAC" ^ AeB ^ ^ 


Kx := 01 


+ 


02 JJ 6(Ao - /i) 
c(Ao - A) 

,,u 

c(A - A„) 


n 


a(Ao — A) 


n 


a(Ao — A) 


h(Xo-X) J-1 A) 

.AeAi.'i ^ ^ AeB ^ ^ ' 


na(A-/z) n 


AeA 


a(A - A) 
t KA - A) 


^n^A-A) n 


'6(A-Ao) 




AgA' 


a(A - A) 
KA - A) 


(2.3) 
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We refer the reader to 23 for a detailed proof. 


Some comments on Theorem]^ are important at this point. For instance, Eq. (2.2) 
is obtained as a particular linear combination of two equations describing more general 
six-vertex model scalar products, i.e. off-shell scalar products. In that case the scalar 
product depends on two sets of var iabl es, i.e. 5°® s^eii ^ ^ofr ^ A„|Af,..., A^), 

and the on-shell case described by (2.2) results from the restriction to fixed Af G B. 


Interestingly, Eq. (2.2) exhibits the same linear structure as the functional equations 


satisfied by partition functions of six-vertex models with domain-wall boundaries 


24,25 . The explicit form of the coefficients (2.3) are the main difference between (2.2) 


and the equations derived in 39,24,25 . In this way one can expect that some features 


discussed in 24 will also hold for (2.2). Before investigating this possibility we present 


some further properties of on-shell scalar products that will be relevant. 


2.1.1 Operatorial formulation 

While on-shell scalar products are functions in n variables, we can readily see that 


Eq. (2.2) runs over the set containing n -|- 1 variables. Thus we have one extra 


variable, namely Aq, which will play a fundamental role for extracting a hierarchy of 


PDEs from (2.2). In fact, this differentiated role played by Aq can be exploited to 


achieve an operatorial formulation of (2.2) with the help of the following definition. 


Definition 4. Write C[Ai’"^, A^^,..., A^^] for the set of meromorphic functions in n 
variables. For each A G A^’"' and any A ^ A^’” define the operator 


Dx ■ C[Af \ ..., A^\..., A^^] ^ C[Af \ ..., A^\..., A^^] 
acting on f G C[Af Af ..., A^^] as 

/) (Ai,..., A,..., An) := /(Ai,..., A,..., An) . 


(2.4) 


(2.5) 


This operator was firstly introduced in 19 and it hts naturally in the structure of 


(2.2). Using this operator we can rewrite our functional equation as £(Ao) 5n(A^’") = 0 
where £ is an operator given in terms of . More precisely, this operator reads 


i;(A„) := A'o + A', A>5" 


( 2 . 6 ) 


AeAi-" 


The importance of this reformulation goes beyond just an optimized notation. For 


instance, this operatorial prescription localizes the whole dependence of (2.2) on Aq into 


the operator £. This feature does not depend on particular realizations of the operator 
which can suitably chosen according to the function space where the desired scalar 
product sits. Interestingly, the operator admits a differential realization with only 
finitely many derivatives in the ring of multivariate polynomials |^. This feature will be 
a key ingredient for our differential description of on-shell scalar products to be discussed 
in Section |2l2] and Section! 
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2.1.2 Properties of Sn 


Theorem [T] has its roots in the Yang-Baxter algebra and carries very little information on 
the particular representation of Bethe vectors we employ to define the scalar product Sn- 
Therefore, one can expect that there will exist different classes of solutions corresponding 
to different representations of Bethe vectors. Linearity is one of the main properties 


of (2.2) and due to that combinations of solutions are also solutions. Thus we need to 


provide conditions capable of singling out the solution corresponding to the desired scalar 
product. This issue is possibly related to the fact that Theorem is not restricted to a 
particular representation of Bethe vectors entering the definition of the scalar product 
Sn- In addition linearity also implies that Eq. (2.2) is at most able to determine the 


desired scalar product up to an overall multiplicative factor independent of the variables 
in the set A^’”. 

Here we are interested in on-shell scalar products of the Wg[s[(2)] six-vertex model 
with twisted boundary conditions 38 , and this choice fixes the representation of Bethe 
vectors entering the construction of Sn- The scalar product Sn for this particular model 
exhibits two important properties that will be required in our analysis. These properties 
are formulated in the following lemmas. 

Lemma 1 (Symmetric function). For 1 < i, j < n let Sij he the permutation operator 
acting on f e C[Af\ ..., by 


(• • • ) . . . , Aj, . . . ) . /(..., Aj, ... , Aj, ...) . 


(2.7) 


Then SijSn = Sn- 


Remark 2. We remark that the symmetry property described in Lemma [I] has been 
already taken into account when writing the functional relation (2.2). It can he shown 
that, reversely, any analytic solution of (2.2) has this symmetry property. 

Definition 5. Let us introduce variables Xi := and set X®’-^ := € C | i < fc < j} 

and := X^’^\{x}. 

Lemma 2 (Polynomial structure). The function Sn is of the form 

SniX^’^) 


SniA^n = 


n 

cexc 


(^- 1 ) 
X 2 


( 2 . 8 ) 


where Sn is a polynomial of degree L — 1 in each of its variables separately. 

The proofs of Lemmas and can be found in 23 . For completeness it is worth 
remarking that, together with the asymptotic behavior of Sn at infinity fixing the leading- 
term coefficient, the above properties uniquely single out the on-shell scalar product 
among the solutions of (2.2). 
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2.2 Partial differential equations 


The possibility of expressing (2.2) as £(Ao)iSn(A^’”') = 0 with operator £ given by (2.6) 


is very snggestive. In fact, we have deliberately chosen to write formula (2.6) as such in 
order to make explicit the resemblance of with some sort of derivative. As we shall 
see, this is not out of reality when the action of is restricted to a particular function 
space. 

Proposition 1. Let C[X^’”] denote the ring of polynomials in n variables xi,X 2 ,... ,Xn 
with complex coefficients. Also, let Cm[X^’"'] C C[X^’"'] be the subspace consisting of 
polynomials of degree at most m in each variable Xj. Then 


= X 


(xo — XiY 


k\ 


dx^ 


(2.9) 


is a realization of (2.5) on the space Cm[X 


’"l.nl 


Proof. The proof follows readily from the Taylor series expansion of arbitrary functions 
in Cm[X^’"']. More details can be found in 19 ,^ 


□ 


Although the realization (2.9) exhibits an appealing structure, the solutions we 
are interested do not immediately belong to Cm[X^’"'] as described in LemmaIn fact, 
the desired on-shell scalar product belongs to the function space := x^“Cl_i[X^’"] 
where x := flsexi." Therefore we can see that the space is not too far from 
Cl_i[X^’"']. On the other hand. Lemmaalso tells us that the function Sn actually lives 
in Cl_i[X^’”']. In this way it is convenient to rewrite (2.2) in terms of functions Sn for 


which the realization (2.9) is available. This motivates the following extra definitions. 


Definition 6 (Polynomial Bethe roots and inhomogeneity sets). Write yi := 


and g := e'*'. In addition to that, it is also convenient to define the sets U := {yi G 


.,2Af 


1 < z < L} and B ;= |xf | 1 < i < n} such that 


n 

yOi 


{xfq - yq-^) (p2 A ^ 


(xf - y) 


n 


xfq ^) 


^xfq 

j¥=i 


Bn-l — xfq) 


( 2 . 10 ) 


Taking into account the above discussion we then rewrite Eq. (2.2) as 

£(i„) S„(X‘’”) = 0 


( 2 , 11 ) 


with operator £ dehned as 


£(xo) := A + K, Dff 
xexi’" 


( 2 . 12 ) 
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The coefficients Kq and in (|2.12) explicitly read 

{xq - XoQ-^) 


Kq ■= (j^iYlixoq- 


yq 




y&A 


n 

.xSX^’" 


{x - Xo) 


n 

xGB 


{xq-xpq 
(x - Xo) 


Kx ■= {q-q ) 


</>2 

y&i 


n 

.xexi'’ 


(xog — xq 
(xo - x) 


n 

x&B 


{xoq — xq 
(xo - x) 


{xo 


X 


X 


h W{xq - yq 

y&A 




n 


xex^ 


{xq — xq 
(x — x) 


-hYiix-y) n 

y&U 


{xq — xq 


X — X 


xGXx 


(2.13) 


The ntility of this reformnlation lies in the fact that £, as dehned by (2.12), now acts 
on Cl_i[X^’"'] and we can readily employ the differential realization (2.9). Thns £ can 
be regarded as a differential operator and its properties will be discnssed in what follows. 

Theorem 2. The operator £ is of the form 


£(xo) = 


{q-q ^)xo 


-I'lT L+n-2 

)Xo 

"•r. - 


Ux&Bi^o - x) 


Xa SI. 


(2.14) 


fc =0 


where VLk are differential operators that do not depend on xq. 


Proof. We hrstly snbstitnte the realization (2.9), adjnsted to the present case, in (2.12) 
and notice that the coefficients Kq and Kx only exhibits simple poles. These poles are 
located at Xq = x* and Xq = xf. It is not hard to see that the residnes of £(xo) at 
the points Xq = Xj vanish identically and, therefore, those points consist of removable 
singnlarities. On the other hand, the residnes of £(xo) at the points Xq = xf only 
vanish up to the Bethe ansatz condition (2.10). Thus we can conclude that £(xo) oc 
11x60(^0 — x)~^ and the remaining part of (2.14) is obtained from the explicit expansion 
of the coefficients Kq and K^ defined in (2.13). □ 


The whole dependence of the equation (2.11) with the variable xq is contained in 
the operator £. As previously remarked, this property is a key ingredient to uncover a 
hierarchy of PDEs governing the scalar product Sn- 

Corollary 1. The set {12^ \ t)<k<L + n — 2] forms a family of PDEs simultaneously 
integrated by Sn: 

flfc = 0 0<A;<L + n-2. (2.15) 

Proof. Consider Eq. (2.11) taking into account Theorem]^ Since Sn does not depend 
on Xo, we are thus left with (2.15). □ 
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Remark 3. The direct inspection of (2.15) for small values of L and n indicates that any 


single equation might already able to determine the desired on-shell polynomial solution up 
to an overall constant factor. This suggests that the set {flk}k might form a commutative 
family of operators, i.e. = 0. Here we shall not investigate this possibility but it 

certainly deserves further attention. 

Although each differential operator Qk can be written down for given values of L and 
n, their explicit form can be rather cumbersome. Fortunately, the leading-term operator 
QL+n -2 exhibits an interesting structure and it can be compactly expressed for arbitrary 
values of L and n. It turns out that the operator flL^n -2 is explicitly given by 


TlL+n-2 = "" - 02 ^” ~ 


X + 


E 


qL 


-1 


yx&B xeXi- 


(L —1)! ^ 

^ ^ xexC" 


(2.16) 


with coefficients 

yx ■= 01 

y&U 




n 


{xq — xq 


[X — X) 


iex^ 


^2\{{x-y) JJ 

y&( xexi’’ 


{xq — xq 
{x — x) 


(2.17) 


The structure of as given by (2.16) and (2.17) exhibits interesting features on 


its own. Firstly, it is similar to the ones found in [^,27,25 for six-vertex models with 
domain-wall boundaries. Moreover, analogously to the PDFs found in 24,27,25 , the 


scalar differential equation = 0 can also be rewritten as a vector-valued first 

order PDF. We shall not perform this construction here as it can be straightforwardly 


obtained from the analysis presented in 27,25 
will be discussed in what follows. 


Solutions of (2.16) for small values of n 


2.2.1 Solutions and their properties 

The inspection of the equation VtLj^ri -2 = 0 for small values of L and n reveals 

that solutions Sn G Ci_i[X^’"'] only exist for parameters xf constrained by Bethe ansatz 


equations (2.10). Our investigations seem to suggest that this single differential equation 


already determines the on-shell scalar product for all L and n. If this indeed the case 


we can readily see from (2.16) that the dependence of 5„(X^’"') with the parameters 
xf should only appear through the combination b := J2xeB^- point out that this 
property is not apparent from Slavnov’s determinant representation of on-shell scalar 
products [^. Next we shall discuss the integration of our PDF for some particular 
values of n. 


Case n = 1. In that case let us first consider L = 2 for which a general analysis 


of (2.16) is within reach. The operator (2.16) then produces the following ordinary 


differential equation. 


i=i i=i 


d5i 

dxi 


+ {(fiq-(j) 2 ){x^ - Xi)Si = 0 , (2.18) 
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whose general solution is given by 


‘5i(a;i) = C exp - log (yl(a;i)) H -arctanh ( ^ 

2 K V K 


(2.19) 


In (2.19) C stands for an integration constant while 


K := ^20102 [q^ivi + y2Y - 2|/i|/2(1 + q^)] - + (t)l){yi - 2 / 2 )^ • (2.20) 

The functions A and u are in their turn defined as 

2 2 

A{x) := 01 JJ(g^x - Vj) - 02?^ 

j=i i=i 

u(x') := ig [ 2 a ;(02 - 0 ig^) + (01 - 02 )(l/i+ 1 / 2 )] • 

Now one can verify that the condition 


uix 


K, 


= 1 


( 2 . 21 ) 


( 2 . 22 ) 


corresponds to the Bethe ansatz equation (2.10). Moreover, assuming (2.22) and recalling 


the identity exp [arctanh (x)] = (1 + x) 2 (l — x)“ 2 , one can check that the solution (2.19) 
collapses to the space Ci[xi]. 

Next we reverse our argument and consider the case n = 1 and arbitrary L starting 
from the assumption 5i G Ci,_i[xi]. In that case is a constant and the integration 


of (2.16) is straightforward. We are thus left with 


5i oc 


~y^ “ y) 

y&U yeU 


(Xi - xf) 


(2.23) 


However, formula (2.23) exhibits a simple pole when xi —>■ xf. Thus, in order to have a 


polynomial, the residue of (2.23) at this pole should vanish. This requirement yields the 
following constraint on xf, 

n =f. ( 2 - 24 ) 


yeu 


(w - y) 


which can be immediately recognized as the Bethe ansatz equation (2.10) for n = 1. 


This analysis suggests that polynomial solutions and Bethe ansatz equations are 
equivalent requirements within our approach. 


Case n 

form 


2. The equation produced by (|2.16|) for n = 2 and arbitrary L assumes the 

(2.25) 


- „n-i g^L-i 


0xf' 
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with coefficients explicitly defined as 


V ■= ^ - (/>2q)(xi + X2 - b) 


Vi ■■= 


V2 := 


(L-1)! 


(L-1)! 


(x2g - Xiq 

(X 2 - Xi) 

(xig - X2g'^ 

(xi - X2) 


y&U 

W{x2g 

y&U 


yg 




yg 




{xiq-X2g 

(xi - X2) 

{X2q - Xiq-^) 

(X2 - Xi) 


n(^> 

y&U 

11(^2 

y&U 


y) 


y) 


(2.26) 


The requirement that 52 ( 2 ^ 1 , 3 ^ 2 ) is a symmetric polynomial living in CL_i[xi,a; 2 ] implies 
a reduction of variables, 


^S2{Xi,X2) 


dxi 


L-l 


= H{x 2 ) and 


^ 52 ( 3 : 1 , 3 : 2 ) 

dx^-^ 


= H{xi) , 


(2.27) 


where now H{x) G Ci_i[a;]. Hence the resolution of (2.25) for S 2 yields the formula 

0~^ [K{xi, X2 )H{xi) - K{x2, Xi)H{x2)] 


S2{Xi,X2) = 

(3:1 - X2) 

with O = (L — 1)! {((>iq^~^ — 02 ?) and 

K{xi,X2) := 0i(3;ig - X2q'^) ]^(3;2g 

y&U 


{xi + X2 - b) 


yg ^) + (t) 2 {x 2 g - Xiq ^)]^( 3:2 

y&A 


(2.28) 


y) • (2.29) 


Interestingly, expression (2.28) can be regarded as a sort of separation of variables in¬ 
duced by the requirement 52 ( 3 : 1 ,X 2 ) G CL_i[a:i,X 2 ]. On the other hand, (2.28) appears 
to have singularities and, by symmetry, it suffices to focus on the variable xi to study 
this issue. Keeping in mind that H{x) is a polynomial, we can see that the RHS of 
(2.28) only exhibits simple poles when xi —>■ X 2 and Xi —)• b — X 2 - As far as the pole 
Xi = X 2 is concerned, the residue of (2.28) vanishes for arbitrary polynomials 77, so it is 
removable. By way of contrast, the residue of (2.28) at the pole Xi = b — X 2 does not 
vanish identically. Thus, in order to have 52 ( 3 : 1 , 0 : 2 ) G £-l-i[xi,X 2 \ we need to require 
the corresponding residue to vanish. This condition implies the following constraint on 
the function 77, 

77(b — x, a:)77(b — x) — K{x,b — x)H{x) = 0 . (2.30) 

Now, since H{x) G Cl_i[x], we can write 77(x) = (2.30) fix the 

coefficients c*. By doing so we immediately notice that this is only possible for values 
of b constrained by a certain polynomial equation. The explicit form of this polynomial 
constraint for L = 2,3 is given in what follows. 


Example 1. For L = 2 we find the constraint 

2 

[g^{fii - 02)b + {yi + y2){hg^ - fii)] JJ ^ ’ 

i=i 


(2.31) 
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while for L = 3 it is given by 


i=i 


(2.32) 


In (2.32) we have, 


T(b) := g^(02 - + 2g^(02 - 0ig)^(0i - 02?) + (02 - 0i?)lVb + Wq , 

i=i 

with parameters W and Wq defined as 


— 0102 


5^2/| + 3 ViVj 

Li=i i<*<i<3 


(1 + Q*) y^yt + y^ 

l<*<i<3 \i=l 


W^o := (01 - 02?') 


f 

3 3 

0102(? + ? ^) 

(?^ - 4 + q~^)yiy2y3 

1 

-1 

1 


+ (01 + 02 ) n yt^ 

l<*<i<3 


(2.33) 


It is interesting to notice that for both examples we hnd that the polynomial determin¬ 
ing b factorizes as P(b) Y\^=i [?^^ ~ + ?^)] ~ 0 where P is a non-trivial polynomial. 

The trivial part consists of a prodnct of linear eqnations, where each eqnation depends 


on a single inhomogeneity i/,. Although those trivial terms are present in (2.31) and 


(2.32), they might not correspond to scalar products of actual Bethe vectors. On the 


other hand, the polynomial equation P(b) = 0 seems to be the relevant one. In fact, one 
can notice its solutions yields the correct number of Bethe eigenvectors. For instance, for 
n = 2 and L = 2 we have only one solution in accordance with the number of admissible 
Bethe vectors. Similarly, for n = 2 and L = 3 we have three solutions which is also the 
correct number of eigenvectors for the six-vertex model in that sector. 


3 Open boundary conditions 


Vertex models of statistical mechanics can be dehned with a variety of boundary condi¬ 
tions and a prominent role is played by lattice systems with open (reflecting) boundary 
conditions. The extension of the QISM to this class of models was put forward in 40 


wherein the so called reflection algebra emerges as a fundamental structure governing 
integrability at the boundaries. Prior to that, the exact diagonalization of certain spin- 


chain hamiltonians with open boundaries had been obtained in 41 42,29 by means of the 


13 

















celebrated Bethe ansatz. In fact, the case of open boundaries required a generalization 
of the Bethe ansatz, which turns out to preserve the scattering theory interpretation. 

The inclusion of this type of boundary conditions within the framework of the QISM 
has several important consequences. For instance, this type of boundary conditions 
have been shown to render systems fully invariant under the action of the underlying 
quantum affine algebra 43 . It is worth remarking here that the study of symmetries 


underlying a physical system is an important step of its analysis. As far as systems with 
open boundaries are concerned, the characterization of the model’s symmetries has even 
found important applications within the context of the AdS/CFT duality 44,45 . In 


addition to that, the extension of the QISM to systems with open boundaries paved the 
way for using a powerful algebraic machinery for studying six-vertex model correlation 


functions 46,47 


The evaluation of correlation functions within the framework of the QISM is inti¬ 
mately associated with the study of Bethe vectors scalar products. More precisely, the 


latter are building blocks of form factors 48,46,47 and this section is devoted to the 


study of scalar products of the six-vertex model with open boundary conditions as so¬ 
lutions of certain PDFs. On the other hand it is worth remarking here that there exist 
alternative methods to study correlation functions for inhnite lattices. In particular, this 


m 


study for open chains was presented in 49,50 using the vertex operator approach and 


5l|"53 by means of the g-Onsager algebra. 


Here our analysis will be based on results recently presented in 26 obtained via the 
algebraic-functional framework. As previously mentioned in Section [23 here we are 
using the terminology on-shell scalar product to refer to the case when only one Bethe 
vector is on shell. That is, the variables parameterizing the Bethe vector are evaluated 
on solutions of the model’s Bethe ansatz equations. On-shell scalar products of the 


six-vertex model with open boundary conditions have been previously studied in 54,46 


using the method of 48 . In that case, a determinant representation resembling Slavnov’s 


formula for the periodic case 16 has been obtained. Although an explicit formula for 


such a quantity is already available in the literature, our interest here lies in possible 
connections with the theory of PDFs and quantum many-body systems. In fact, similar 
connections have been observed for several quantities tackled through the algebraic- 
functional method. For instance, this is the case for partition functions of six-vertex 
models with domain-wall boundary conditions 19,24,25 and the spectral problem of 
vertex models transfer matrices (271. In what follows we shall demonstrate that this 


is also the case for on-shell scalar products associated with the six-vertex model with 
diagonal open boundary conditions. 


3.1 Functional equation 

Among the main results of we have a set of two functional equations characterizing 
scalar products of Bethe vectors for the six-vertex model with open boundary conditions. 
The simultaneous resolution of those functional equations, taking into account certain 
properties expected for off-shell scalar products, yields the desired solution in terms of a 
multiple contour integral. Similarly to the case with boundary twists considered in 23 


the on-shell case for open boundaries can also be described by a single equation obtained 
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as a linear combination of the equations presented in j^. Here we shall make use of 
dehnitions and conventions introduced in since our present analysis builds on results 
of the aforementioned work. Moreover, the open boundaries case will also require the 
following notation, in addition to definitions [T] and given in Section 

Definition 7 (Double-row Bethe roots). Define the set S := {Af, Af,..., A®} for e C 
solving the following system of Bethe ansatz eguations for 1 < k < n, 


6(Af + h) 6(Af - h) 
a(Af - h) a(Af + h) 


n 




a(Af -/i) a(Af +/i) 

- h) + h) 


ir'n 


®(-^f ~ -^f) ®(-^f + Xf + j) 


1= 

i^k 


, (i(Af - Af) b{\f + Af) 


(3.1) 


Theorem 3 (Functional equation). Let Tn- t C denote on-shell scalar products 
for the six-vertex model with open boundary conditions. Then Tn satisfies the functional 
eguation 

Lo TniXfin + r„(A°/) = 0 , ( 3 . 2 ) 

AeAi." 

with coefficients Lq and Lx defined as 

a(2Ao 7 ) 


Lq 6(Ao -l- h)6(Ao — h) 


X 


n 


ci(A — Aq) &(A -|- Aq) 


a(Ao — p.)a(Ao -I- /i) 

ci(A — Aq) &(A -|- Ao) 


n 


LAeAi'" ~ Aei ~ 


-|- ci(Ao — h)a(Ao -l- h) 

n 


6 (2Ao 


X 


.AeAC" 


®(Ao — A) ci(Ao -l- A -|- 7 ) 
b{Xo - A) 6 (Ao + XX-'j) 


JJ 6(Ao -/i)&(Ao-h/i) 

a(Ao - A) a(Ao -h A 7 ) 


n 

xes 


b{Xo - A) b{Xo A 7 ) 


(3.3) 


Lx := 


a(2Ao -l- 7 ) c(Ao — A) 6(2A) 
a(Ao -l- A) 6(Ao — A) a(2A) 


X 


b{X h)b{X — h) a(A — p)a{X -|- /i) 




AgA: 


a(A - A) hiX + A) 
„ h(X — A) a(A -I- A) 


i(A — h)a{X -\- h) 6 (A — /i) 6 (A -I- p) 


f-L^lA 


AgA 


(z(A — A) ci(A -|- A -|- 'y) 
KX - A) 6 (A + A + 7 ) 


Proof. Considering the results of 26 , we multiply Eq. type A by b{h — Aq) ^[ 2 ^ 0 + 7 ) 
add it to Eq. type D multiplied by a(Ao -l- h). Then we are left with (3.2) by assuming 
that the variables Af are constrained by the Bethe ansatz equations (3.1). □ 


15 


























tional relations derived in 39 


It is interesting to notice that (3.2) exhibits the same structnre as (2.2) and the fnnc- 


24,25 for partition fnnctions with domain-wall bonndaries. 


In this sense, both (2.2) and (3.2) seem to be part of a general structure accommodating 


several quantities associated with integrable vertex models. Next we shall discuss the 


passage from (3.2) to a hierarchy of PDEs. 


3.1.1 Operatorial formulation 


The procedure described here will go along the lines of that presented in Section |2.1.1 
We shall start by reformulat ing t he functional equation described in T 
of operators D^° dehned in (2.5). This is the hrst step for converting i 


leor em in terms 


3.2) into a family 


of PDEs. In fact, this reformulation only depends on the structure of the equation and 


we can readily notice that (3.2) exhibits the same structure as (2.2). Therefore, we shall 


restrict ourselves to presenting only the hnal results. 
Let OH be an operator dehned as 


'm(A„):=L„+ L, 




(3.4) 


AeAi>" 


In terms of fOT the functional equation (3.2) can be written as IDl(Ao) 7^(A^’ ^) = 0. 
One important aspect of this reformulation is that the operators entering ( 3.4[ ) are 
dehned in a rather abstract way (2.5). This feature leaves room for choosing suitable 


realizations depending on the function space accommodating the desired on-shell scalar 


product Tn- This will be one of the aspects discussed in Section 3.1.2 


3.1.2 Properties of 7^ 


Most of the discussion of Section 2.1.2 applies to Eq. (3.2) and here we restrict our¬ 


selves to presenting only special properties expected from the scalar product Tn- Those 


properties are meant to guide us through the resolution of (3.2) ensuring that we have 


actually selected the solution corresponding to on-shell scalar products associated with 

In fact, the properties we shall 


26 


the Uq[s{{2)] six-vertex model with open boundaries 
use to single out actual on-shell scalar products characterize a function space; and for 
such function space ( 3.2[ ) gives rise to a hierarchy of PDEs. The required properties will 
be stated in lemmas ^ andwhose proofs are available in [2^ . 

Lemma 3 (Symmetric function). The function T is symmetric under the permutation 
of variables, i.e. SijT = T for all I < i, j < n. 

Remark 4. As for the case with twisted boundary conditions, this symmetry property is 


already manifested in the functional relation (3.2) 


Lemma 4 (Polynomial structure). The function T is of the form 

7;(xi'-) 


Tn{A^n = 


n 

xexi’" 


(3.5) 


where T G C 2 l[X^’"'] is a multivariate polynomial of degree 2L in each variable separately. 
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3.2 Partial differential equations 

Solutions of 9Jt(Ao) 7^(A^’"') = 0 corresponding to on-shell scalar products are character¬ 
ized by Lemmas and 1^ In other words, the desired scalar products consist of solutions 
of our equation in the space Qn ■= x“^C 2 i,[X^’"']. The situation here is analogous to the 
case discussed in Section 2^ and a differential description of on-shell scalar products 
for the case with open boundaries is available in terms of functions Tn- This is due to 
the fact that Tn G C 2 l[X^’"'] for which (2.9) can be employed. The following additional 
dehnitions will also be useful. 

Definition 8 (Polynomial double-row Bethe roots). Considering Definition^we further 
write t := and t := e^. We then define 8 : = \x f, xf ,... ,x^} formed by variables 
xf = e'^^k G C satisfying the polynomial form of (3.1) for 1 <k <n, namely 


(xft - t 1 ) (xft ^ -t) -p-p (xfg - yq i) (xfg - y T 


n 


{xft-^q - fg-i) {xftq - W^q-^) (xf - y) 


(4 - y4 


n 


xf g - xfq (xf g2 ^xf) ^q ‘^) 


(3.6) 




xfq) 


(xt 


(4)4 


Next we rewrite (3.2) in terms of functions T described in Lemma With the help 
of Dehnition 1^ we are then left with the equation IH(xo) 7)i(X^’"') = 0 with operator Tft 
dehned as 

im(xo) := To + 5^ 4 DT . (3.7) 


The coefficients Lq and in (3.7) are in their turn given by 


Lo := (4^ -XqC ^)(x4 


,2+-l 


Yl(xoq-yq 4xoq-y ^) 


X 


n 


(xog-xo^g-^) 
(xg - xog"^) {x - Xq j-j 


y&A 

(xg - xog"^) ( 


X — Xr 


4 ) 


(x-xo) (a;g-Xo^g-^) 
(4 q4 - x4q-h) (x| qt - x4q-H4 


x£E 

{xo 


(x - Xo) (xg - Xq ^g-i) 


X 


(xog-xo^g-i) 


r4lY “ 


y 


-1^ 


X 


n 

.a:eXl’' 


(xog — xg ^) (xog^ — x ^g ^) 
(xo — x) (xog —x“^g“^) 


n 

x££ 


y&A 
^-1 


{xoq — xq ^) (xog^ — X ^g ^) 
(xo — x) (xog — x“^g“^) 


(3.8) 
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Lx • ^0 


{q-q {xoq^ - Xp^q (x-x 
(xo — x) (xoq — x~^g~^) (xg — x~^g~^) 


X 


(x^t — x~ 2 t~^)(x^t~^ — x~^t) ]^(a:g — yg~^){xg — y~^g~^) 

y&A 

{xg — xg~^) {x — x“^) 


X 


n 


xexi 


{x — x) {xg — x~^g~^) 


{x^gt ^ — x 2g ^t){x‘2gt — x ^g ^)'^^{x — y){x — y 

y&U 


X 


n 


xexi 


{xg — xg (xg^ — x ^g 
{x — x) {xg — x~^g~^) 


The operator VJt now acts on C 2 l[X^’”] and we can simply use the differential realiza¬ 
tion (2.5) adapted to this case. Thus can also be regarded as a differential operator. 


Its structure is described in the following theorem. 
Theorem 4. The operator 50T is of the form 


DJl{xo) = 


(q-q 


2L+3n 


JJ (xog-x ^g ^)JJ(xog-x ^g ^){xo-x) fr 




(3.9) 


a:eXC" 


x&£ 


In (3.9) we use <I>fc to denote differential operators independent of Xq. 


Proof. We proceed by inserting the differential realization (2.5) into (3.7). Then we notice 
from (3.8) that fOT, as a function of Xq, only contains simple poles. These are located at 
Xo = xo G X^’", Xq ^ G g^ X^’", xq G ^ and Xq ^ G E. The residues at the poles 

Xo = ±g“^ and xq G X^’” vanish identically, hence they consist of removable singularities. 
On the other hand, the residues at the remaining poles do not vanish identically and 


give rise to the denominator in (3.9). The numerator is now a polynomial expression in 
Xq. The degree of this polynomial can be determined as follows. In order to compensate 


the dependence with xq in the denominator of (3.9) we need to include a factor of x; 


^3n 


in the polynomial expansion. Next we count the powers coming from and notice 

that the leading term is proportional to Xq^. Although the leading term coming from Lq 
seems to contain one more power of xq, its careful inspection shows that the coefficient 
of that term actually vanishes. Thus Lq also yields a leading term proportional to Xq^. 
Thus we are left with a polynomial in xq of degree 2L -|- 3n. □ 


Our analysis will now proceed along the lines described in Section 2^ For that the 
reformulation of (3.2) as IH(xo) 77(X^’"') = 0 plays a crucial role. 


Corollary 2. The set of differential operators {4*^} constitutes a hierarchy of PDFs 
satisfying 

4)fc 7;(X^’^) = 0 for 0<k<2L + 3n. (3.10) 
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Proof. The proof is the same as for Corollary □ 

Similarly to the case with twisted boundaries, the leading-term coefficient ^ 2 L+zn 


in the expansion (3.9) also plays a special role. This coefficient can be written down 
explicitly for arbitrary values of n and L in a compact form. More precisely, it reads 


$ 


2L+3n 




c — {xq + X ^q 

xSXl-" 


+ 


Q 


271+1 


^ {x — X ^ 8“^^ 


(2L)! {xq — x~^q~^) ^dx'^^ ' 


(3.11) 


where the parameter c ;= ^ contains the whole dependence of ^ 2 L+ 3 n 

with Bethe roots. In its turn the function has been dehned as 

■= {x'^t — x~' 2 t~^){x 2 f~^ — x~^f)Y\_{xq — yq~^){xq — y~^q~^) 

y&A 

{xq — xq~^) {x — x~^) 


X 


n 




{x — x) {xq — X ^q 
— {x^qt~^ — x~^q~^t){x^qf — x~^q~^t~^) — y){x — y~^) 

{xq — xq~^) {xq^ — x~^q~‘^) 


y&A 


X 


n 


*exi 


{x — x) {xq — x~^q~^) 


(3.12) 


From (3.11) we can readily see that equation ^ 2 L+ 3 n = 0 is a linear PDF 


of order 2L in n variables. Its structure is quite similar to (2.16) and the PDFs found 


m 


Z7\ 25 describing integrable vertex models. Although the order of (3.11) depends on 


the lattice length L, the corresponding differential equation can also be recasted as a 


hrst-order vector-valued PDF along the lines of 27 . Moreover, the fact that Bethe roots 


only enter (3.11) through the parameter c will play an important role for its resolution. 


Section 3.2.1 will then be devoted to the analysis of the equation ^ 2 L+ 3 n = 0 


for small values of the magnon number n. 

3.2.1 Solutions and their properties 

Similarly to the case with boundary twists discussed in Section a determinant repre¬ 


sentation for Tn also exists 54,46 . However, it is again unclear from those results that 


Tn might depend on Bethe roots only through a particular combination. In this section 
we shall focus on the resolution of the single equation ^ 2 L+ 3 n = 0, although 

Corollary entitles us with 2L -|- 3n -|- 1 PDFs. A priori, it is not clear that a single 


equation would be able to £x the desired solution but direct inspection of (3.10) for small 


values of L and n suggests that each equation separately is able to determine the de¬ 
sired on-shell polynomial solution up to an overall multiplicative factor. If this is indeed 
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the case then our results imply that the on-shell scalar products only depend on the 
Bethe roots through this particular combination c. Moreover, we also hnd that solutions 
Tn G C 2 l[X^’"'] only exists for parameters xf constrained by the system of equations 
(3.6). In what follows we shall demonstrate this fact explicitly for the cases n = 1, 2. 


Case n = 1. In this case our problem consists of solving the following ordinary differ¬ 
ential equation, 


(xi - ^ d^^Ti 

(2L)! (xig — dxf^ 




t ^tq^) (c - Xiq - Xi 7i = 0 . 

(3.13) 


For solutions 71 G C 2 l[xi] we have that *^^ 2 ? is a constant and we are left with 


Ti oc 


(xi - Xi 




(xig 


Xi ^q (xig 


Xi ^q ^ 


(3.14) 


Formula (3.12) for n = 1 shows that contains no poles for finite values of the variable 
Xi. Hence, the only source of poles in (3.14) is the denominator of (3.13). More precisely, 
these poles are located at 

i. (xig)^ = 1 , 

a (xig)^ — Xiq c -|- 1 = 0 . 


The condition 71 G C 2 l[xi] now asks for the residues of (3.14) at those poles to vanish. 
As far as the poles (z) are concerned, we can readily see that the corresponding residues 
vanish identically. However, the analysis for the poles at {ii) is slightly more involved. 
For that we hrst recall the dehnition of c and notice that {ii) can be rewritten as (xi — 
xf)(xig^ — 1/xf) = 0. The residue of (3.14) at xi = xf vanishes provided Z^b = 0. This 
condition is equivalent to the Bethe ansatz equation (3.6) for n = 1. Likewise, the pole 
at Xi = l/(g^xf) has zero residue when xf is on shell. The latter follows from the fact 
that the transformation xf i—)■ l/(xfg^) only modihes (3.6) by an overall factor. Hence, 


our analysis so far suggests that polynomial solutions and Bethe ansatz equations are 
equivalent requirements for the equation 4>2L+3n ’^(X^’"') = 0 also for the case of open 
boundaries. 


Case n = 2. We proceed to the case n = 2 along the same lines described in Sec¬ 
tion 


2.2.1 For that we assume T 2 G C 2 l[xi,X 2 ] and from (3.11) one then finds 

.A(xi,X2).H(xi) - .A(x2,xi).H(x2) 


7^(Xi,X2) = 


0-1 


(xi - X 2 )(xiX 2 g - g i) [(xiX 2 )(xiX 2 g + g i) - xiX 2 c] 


(3.15) 
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where O := (2L)! {tt ^ — t The function H in (3.15) is a polynomial, more 

precisely H{x) G C. 2 l[x\-i whilst K is dehned as 


K{xi,X2) ■■= 


Xi (x2 - X 2 
(x2q-X2^q-^) 


X 


(X2t-t ^)(x2t ^ - t)(xiq - X2q ^){xiX2 - l)W{x2q - yq ^){x2q-y ^q 

y&A 

+ {X2qt~^ - q~^t){x2qi - q~H~^){x2q - Xiq~^){xiX2q^ - q~‘^) 


X - y){x2-y 

y&A 


(3.16) 

In accordance with Lemmaformula (3.15) is manifestly symmetric under the permu¬ 
tation of variables Xi and X 2 . Th is pro perty allows us to focus on the variable Xi in order 


to analyze the poles structure of (3.15|). The function H, although not yet specified, does 
not contribute to this poles struct 


ure as it is a polynomial. Thus the RHS of (3.15) only 


exhibits simple poles at 


i. Xi = X2 , 
a. Xi = -l/{x2q^) , 

in. (xig)^ = 1 , 

iv. (xi -t- X 2 ){xiX 2 q + q~^) — X 1 X 2 c = 0 . 


The residue at each of the poles {i)-{iii) vanishes, so these singularities are removable. 
By way of contrast, the poles at (iv) require a more detailed analysis. In order to carry 
out this analysis we introduce parameterizations u := xi +X 2 and v := 0 : 1 X 2 . By doing so 
we avoid the appearance of square roots when looking at the algebraic curve defined by 
{iv). The requirement that the residue of ( |3.15[ ) vanish at {iv) then yields the constraint 

k{xi{u,v),X2{u,v))H{xi{u,v)) - k[x2{u,v),Xi{u,v))H{x2{u,v)) =0 (3.17) 


for u = vc{vq + q~^)~^. Next we write H{x) = X]i=o (3.17) fix the coefficients 

Cj. At this point it is important to emphasize that the use of variables u and v plays 
a fundamental role. They overcome the problem of dealing with square roots in the 
analysis of (3.17) allowing one to read out a linear system of equations for the coefficients 
Cj. Similar to the case with twisted boundaries discussed in Section 2.2, here our linear 
system also has more equations than coefficients Cj to be fixed. The extra equations 
imposes a constraint on the parameter c that will be explicitly written down for L = 2, 3 
in what follows. 
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Example 2. For L = 2 we find the following condition on the parameter c, 


( - q-H^){yi + ?/2)(l + l/il/2) \ 

V - <it‘^)yiy2 ) 

2 

n + i/r^)) = 0. 

i=i 


(3.18) 


In order to avoid long formulae, for L = 3 we restrict ourselves to the homogeneous case 
yi = 1. In that case we find that c has to satisfy 


(^(t 2 _ ^3 _ ^^2 _ ^2 ^ c - IVo) n “ 2(g + q-^)) 

i=i 


with coefficients Wo, Wi and W 2 explicitly reading 


0 , 

(3.19) 


^2 := {q - q~^){q~^ + A + q^){l + t^P) + 12(g-V - qP) , 

W, := (q - q-^f(2 + q^)(2 + q-^)(l + fH^) 

+ 6(q - q~^) [(2 + q^)(2 + q~^){q~H^ - qP) - + g^P)] (1 + t‘^P) , 

- (g^ - 51g“^ + 2g-=^)g-H^ - {2q^ - 51g“^ + g“^)gt^ 

+ (g® + 4g^ - 13g2 - 80 - 13g-2 + 4g-^ + q-^)t‘^P 

^o:= {q-q-^){q"-q-")\l + t^^) 

+ 6(g - g-^)(g2 - g-^) [(2 + q-^)t^ - (2 + q^)P] (1 + tH^) (3-20) 

- - 37g - 13g-^ + q-^)q-H^{l + 

- {q" - q~^W - 13g - 37g-i + g-^)g¥(l + t^P) 

+ {q" - q~"){q + q~^){q" + g' - 52 + g-^ + g-")t¥(i + 

+ 2(3 + g-2)(g2-12 + 3g-2)g37® + 2{3 + q‘^){3q‘^ - 12 + q-^)q-H^ 

+ 2(16g^ - 36g^ - 57g - 31g"^ + 9g"^ + 3g"^)t^? 

- 2(3g® + 9g^ - 31g - 57g“^ - 36g“^ + . 


Remark 5. The structure of the polynomial conditions (3.18) and (3.19) is essentially 


the same as the one found for the parameter b in Section \2.2. 1\ Thus one can see that the 
counting of admissible solutions c is also in agreement with the number of eigenvectors 
for the six-vertex model with open boundary conditions. 


4 Concluding remarks 


The evaluation of form factors for quantum integrable systems was shown in 16 to 
be directly related to the computation of vertex models scalar products, and in this 
work we have presented a description of such quantities in terms of linear PDEs. It is 
worth remarking that the description of certain correlation functions in terms of linear 
PDEs, namely KZ equations, is among the achievements of conformal symmetry [^. In 
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this sense we can see a convergence of those approaches at least at an ideological level. 
Althongh it is not clear if there exists a precise relation it wonld be interesting to fnrther 
investigate this possibility. 

In the present paper we have selected two types of vertex models to illnstrate onr 
method: the six-vertex model with generalized toroidal bonndary conditions, also known 
as twisted bonndaries; and the six-vertex model with open bonndary conditions. The 
twisted bonndaries case has been discnssed in Section |2] while Section |3] is concerned with 
the open bonndaries case. It is also important to stress here that our differential approach 
is one of the outcomes of the algebraic-functional method for scalar products proposed 
*2^ and (26|. Within this framework, the Yang-Baxter and reflection algebras plays 


m 

a fundamental role, and the PDEs obtained in Section |2.2| and Section |3.2| are a direct 
consequence of those algebraic structures. This feature is also in consonance with the 
derivation of KZ equations by means of vertex algebras. However, our approach not 
only renders a single PDE describing on-shell scalar products, but it produces a whole 
hierarchy of PDEs for each one of the scalar products considered here. 

Our PDEs also unveil interesting properties of on-shell scalar products for six-vertex 
models that were not apparent from the results previously obtained in the literature 
From (2.22) and (3.11) we can see that the dependence of iS„ and Tn with 


the corresponding Bethe roots might only appear through particular linear combinations 
denoted respectively by the parameters b and c in the main text. In Section 2Y we have 
shown that b is determined by a single polynomial equation. This appears to represent a 
signihcant simplihcation since the standard approaches require the resolution of a system 
of Bethe ansatz equations (2.10). In this way, our approach somehow seems to be able to 
combine the whole set of Bethe ansatz equations into a single equation for the quantity 
b, which is the one relevant for the on-shell scalar product Sn- The same applies to the 


parameter c in the case of Tn as shown in Section 3.2 


Nevertheless, there are still many questions that have eluded us so far. Foremost 
it would be important to prove rigorously that the leading term differential operators 
QL+n -2 and ^ 2 L+ 3 n are indeed able to fix the desired scalar products. The particular 
combinations of Bethe roots forming the parameters b and c are quite appealing and one 
might wonder if there is a more concrete meaning associated to them. Those param¬ 
eters consist of simple sums of the associated Bethe roots, which suggests they might 
correspond to the eigenvalue of one of the operators contained in the commuting family 
generated by the vertex model’s transfer matrix. Moreover, it is unclear to us at the 
moment how the system of Bethe ansatz equations could be suitably combined yielding 
the single equation for the relevant parameter b or c. Those questions certainly deserve 
further investigation. 
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